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We propose a generalization of the supersymmetric representation of spins with symplectic sym¬ 
metry, generalizing the rotation group of the spin from SU{2) to SP{N). As a test application of 
this new representation, we consider two toy models involving a competition of the Kondo effect 
and antiferromagnetism: a two-impurity model and a frustrated three-impurity model. Exploring 
an ensemble of L-shaped representations with a hxed number of boxes in their respective Young 
tableaux, we allow the system to choose which representation is energetically more favorable in 
each region in parameter space. We discuss how the features of these preliminary applications can 
generalize to Kondo lattice models. 


I. INTRODUCTION AND MOTIVATION 

Heavy fermion materials involve a lattice of localized 
magnetic moments derived from f-electrons, embedded 
in a conduction sea formed principally from delocalized 
d-electron^^. The physics of these materials can be un¬ 
derstood as a consequence of the interplay between two 
competing physical processes: 

- the Kondo effect, which tends to screen the local 
moments to produce a band of heavy electrons, and 

- antiferromagnetism, which locks the local moments 
together via the RKKY interaction, into a state 
with long range magnetic order. 

The characteristic scales for these two processes are 

Tk ^ Trkky ^ PqJk, (1) 

where Jk is the strength of the onsite Kondo interaction 
between the localized f-electrons and conduction elec¬ 
trons, po ^ 1/7^ is the density of states of the conduction 
electrons at the Fermi energy and D the bandwidth. The 
Kondo temperature, T^, sets the energy scale for the on¬ 
set of the Kondo effect and consequently the formation 
of a coherent heavy Fermi liquid (HFL) while Trkky 
defines the energy scale for the onset of magnetic order. 

The family of heavy fermion materials provides an im¬ 
portant setting for the study of quantum criticalit^EEl 
which develops when a continuous second-order phase 
transition is suppressed to absolute zero temperature. 
The small characteristic energy scales of these com¬ 
pounds makes them highly tunable, allowing the ready 
exploration of the phase diagram as a function of pres¬ 
sure, magnetic field or doping. Superconductivity is of¬ 
ten found in the vicinity of magnetic quantum critical 
points (QCP). At temperatures above the quantum crit¬ 
ical point non Fermi liquid (NFL) behavior is observed, 
generally characterized by sub-quadratic temperature de¬ 
pendence of the resistivity p and a logarithmic tempera¬ 
ture dependence of the specific heat coefficient 7 = ^, 

pocT^, (a<2) 

1 . (2) 


where Tq is the characteristic scale of the spin fluctua¬ 
tions. For a review of experimental properties of these 
materials see StewariP. 

One of the central challenges of heavy fermion materi¬ 
als is to understand the mechanism by which magnetism 
develops within the heavy electron fluid. Traditionally, 
magnetism and heavy fermion behavior have been re¬ 
garded as two mutually exclusive states, separated by 
a single quantum critical point. However, a variety of 
recent experiments suggest a richer state of affairs, in 
particular: 

- YbRh 2 Si 2 can be driven to a quantum critical point 
by the application of magnetic field, where both 
the Neel temperature and the Kondo energy scale 
appear to simultaneously vanish. However, when 
doped, these two energy scales appear to separate 
from one-another, indicating that the break-down 
of Fermi liquid behavior and the development of 
magnetism are not rigidly pinned togetheJ^^® 

- In the 115 superconductor CeRhIns there is evi¬ 
dence for a microscopic and homogeneous coexis¬ 
tence of local moment magnetism and supercon¬ 
ductivity under pressur^ 

- Neutron scattering experiments observe a par¬ 
tially ordered state in the geometrically frustrated 
CePdAl, in which one third of the Ce moments 
do not participate in the long-range order, sug- 
gesting the development of inhomogeneous Kondo 
statej™. 

Various phenomenological frameworks have been pro¬ 
posed for the understanding of heavy fermion sys¬ 
tems. The classical framework proposed in the 70’s 
by DoniaclP^, involves a competition between Tk and 
Trkky determining the ground state to be a heavy Fermi 
liquid or magnetically ordered. More recently a new axis 
was added to this picture, by the inclusio n of g eomet- 
ric frustration or reduction of dimensionalityT^Eil These 
two factors contribute towards the suppression of mag¬ 
netism in a different way, if compared to the competition 
with the Kondo effect. Also, based on experiments in 
several families of heavy fermions, a phenomenological 
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two-fluid picture was proposed by Nakat suji a nd Pines, 
with predictive power on the ground state^^n^. 

Unfortunately these proposals do not give us infor¬ 
mation about the character of the transition between 
the HFL and magnetic phases, and its theoretical de¬ 
scription has remained an unsolved challenge for several 
decades. Theoretical proposals based on a spin density 
wave description of the Kondo breakdowiPIl, 

deconfined quantum criticality^ and local quantum 
criticalit 3 ^^ have been suggested, but no one picture is 
yet able to fully account for experimental observations. 

A. Spin Representations in the Kondo Model 

The Kondo lattice Hamiltonian 

Hkl = ^ ekcj^^Ckcr + Jk ^ Si • Si, (3) 

kcr i 

provides a minimal model for heavy fermion systems. 
The first term in Hkl describes a band of conduction 
electrons with dispersion Ck, Jk is the antiferromagnetic 
Kondo coupling between the local moment Si and the 
spin of the conduction electron Si at site i. 

The local moments are neutral entities uniquely char¬ 
acterized by their spin quantum numbers. The removal 
of the charge degrees of freedom from the Hilbert space 
of the localized f-electrons means that spin operators 
do not follow canonical commutation relations; conse¬ 
quently, their treatment within a path integral or dia¬ 
grammatic approach is complicated by the absence of a 
Wick’s theorem. To circumvent this difficulty, the spin 
operator is traditionally factorized in terms of creation 
and annihilation operators: 

SaH = (4) 

where are bosonic or fermionic creation and an¬ 

nihilation operators, respectively, and the indexes f3 = 
{1,2} for an SU(2) spin. There are actually sev¬ 
eral such spin representations: the Holstein-PrimakoflP^, 
Schwinger bosorP^, Abrikosov pseudo- fermioiP^ and the 
drone or Majorana fermiorP^ representations, among 
others. 

The physics that each of these representations de¬ 
scribes is profoundly different. For example, the anti¬ 
ferromagnetic (AFM) phase at small Jk is very effec¬ 
tively described by a Schwinger boson representation of 
the local moments, with the condensation of t he bo sons 
corresponding to the onset of magnetic ordei^^^^. By 
contrast, the heavy Fermi liquid phase at large Jk is 
successfully captured by a fermionic representation of the 
spinP^. We take the view that the success of these two 
representations in the different limits is not simply one 
of mathematical convenience; rather, it reflects the phys¬ 
ical transformation of both the spin correlations and the 
excitations of the local moments: these evolve from col¬ 
lective spin waves to charged heavy fermions. Remark¬ 
ably, experiment indicates that these two phases connect 


together continuously via a quantum critical point, sug¬ 
gesting that at quantum criticality the two representa¬ 
tions merge. 

In this paper we argue that a full description of heavy 
fermion materials requires a methodology that can cap¬ 
ture the transformation in the character of the ground- 
state and its spin excitations. This, in turn, leads us to 
adopt a supersymmetric representation of the spiiP^ 

Sa0=flf0 + bib0. (5) 

Here, /j^, fa and ba are respectively, fermionic and 
bosonic creation and annihilation operators. The spin 
is supersymmetric because it is invariant under transfor¬ 
mations that take bosons into fermions and vice versa; 
these are generated by fermionic operators which will be 
introduced in the next section. 

One of the challenges of such a factorization, is that it 
requires a constraint which guarantees that the physics 
lies within the physical Hilbert spac^^. For example, an 
elementary spin S' = 1/2 Kramers doublet requires the 
constraint Q = nij ^ Uf = 1. Within this constrained 
Hilbert space, the most general wavefunction is an en¬ 
tangled product 

\^)=Pg{\^f)<S)\^b)), ( 6 ) 

where |T^) and |Ti?) are the bosonic and fermionic com¬ 
ponents of the wavefunction, respectively, while Pq is 
a Gutzwiller projection operator. This operator can be 
written as: 

Pg = /ni —(7) 

J 27r 

which imposes the constraint nsi P npi = 1 at each site 
i. The unprojected wavefunction describes the for¬ 
mation of long-range magnetic correlations in the form 
of a bosonic RVB wavefunction, while \ fjF) captures the 
development of Kondo of singlets and the development of 
a large Fermi surface of heavy electrons. The Gutzwiller 
projection entangles the two components of the wave- 
function into a single entity as illustrated in Fig. 



FIG. 1. Schematic illustration of a Gutzwiller wavefunction 
formed from the projected product of a bosonic RVB wave- 
function and a Kondo-screened Fermi liquid, to form an en¬ 
tangled combination of both wavefunctions. 
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This merged wavefunction has, in principle, the poten¬ 
tial to capture the two-fluid aspect of the heavy fermion 
ground-state. 


B. Large-Approach 

The absence of a small parameter in the Kondo model 
effectively rules out the use of conventional perturba¬ 
tion theory. The alternative approach, followed here, 
is the use of a large-A expansion in which the funda¬ 
mental representations contain A, rather than 2 compo¬ 
nents. In this approach 1/N ^ hs plays the role of syn¬ 
thetic Planck’s constant leading to a controlled mean- 
field (“classical”) theory in the large-A limit, with the 
possibility of expanding the fluctuations and the con¬ 
straint condition as a power-series in 1/N about the 
large-A limit. The simplest generalization takes SU{2) 
to Written in traceless form the SU{N) 

spin is then 

= flu + - {np + ub) Sap/N, (8) 

where a, /3 G {1, 2,..., N}. However, in this paper we seek 
to extend the supersymmetric description of spins to the 
symplectic subgroup SP{N) of SU ( 

SlN'’ = flf^ + ^a^/5 “ + bli^b_a)- (9) 

Here N must be even, while the range of the elemen¬ 
tary spin quantum numbers is a, P G {±l,±2,...,±7V/2}. 
The tilde notation, employed extensively in this article, 
denotes the sign of the index 

a = sgn((a), (10) 

with the analogous definition for other indexes. This 
new spin operator has the symplectic property 5'f = 

—(and is thus also traceless). The symplec¬ 
tic group SP{N) offers many advantages for condensed 
matter physics, allowing for a consistent extension of the 
notion of time reversal symmetry to the large-A limit, 
which permits one to form singlet pair s of p articles that 
are absent in the SU{N) generalizatioiP^EH. This capa¬ 
bility is vital to describe antiferromagnetism and super¬ 
conductivity. 


w 



FIG. 2. Showing an L-shaped Young tableau, characterized 
by two parameters, the width w and the height h of the 
tableau. 



Bosonic 

(magnetic) 


Mixed 


Fermionic 

(paramagnetic) 


FIG. 3. Series of i tableaux in the large-Y limit ranging from a 
fully symmetric (top) towards a fully antisymmetric (bottom), 
passing through a series of L-shape representations. 


The concept of a supersymmetric spin was intro duced 
in previous studies of impurity Kondo mo del In 

the work presented here, we follow the lines of Coleman 
et al.l^, with the additional generalization to 5'P(Y), 
and discuss the spin representations in terms of Young 
tableaux. Young tableaux provide a precise pictorial 
rendition of irreducible spin representations: horizontal 
Young tableaux label completely symmetric representa¬ 
tions, which are naturally described by bosons, while ver¬ 
tical Young tableaux label completely antisymmetric rep¬ 
resentations, usually described by fermions. The use of 
supersymmetric representations lead us to consider the 
set of representations characterized by L-shaped Young 
tableaux (Fig. [^. These representations are character¬ 
ized by two constants: 

- the total number of elementary spins (or boxes) in 
the representation Q = h P w — where h and w 
are height, and width of the Young tableau, respec¬ 
tively, and 

- the asymmetry Y = h — w of the L-shaped Young 
tableau, as discussed in Coleman et al.EH. 

The asymmetry of the representation is absent in a 
physical SU (2) spin-1/2, in which case the Young tableau 
is depicted by a single box, but once we enlarge the sym¬ 
metry group of the spin in order to develop a large-Y the¬ 
ory, we find a family of representations that range from 
a completely symmetric representation, fully described 
by bosons, to a completely antisymmetric representation, 
described only by fermions, including a whole plethora 
of intermediate representations that we refer to as mixed 
representations^ depicted by L-shaped Young tableaux 
(see Fig. [^. The possibility of mean-field solutions de¬ 
scribed by mixed representations is interesting as it may 
permit the description of new states of matter, including 
coexistence of magnetism with superconductivity or with 
heavy Fermi liquid phases. 
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For a given value of one needs to decide which rep¬ 
resentation to choose in order to proceed with the calcu¬ 
lations. Traditionally a purely bosonic representation or 
a purely fermionic representation is chosen, but the su¬ 
persymmetric approach provides the possibility of con¬ 
sidering an L-shaped representation. To constrain the 
problem to such a representation one must fix the val¬ 
ues of Q = Qo and Y = Yq through the introduction of 
projection operators into the partition function: 

Z = Tr[PQ,,Yoe-^^]. ( 11 ) 

Typically, the more negative T, the more symmetric 
the spin representation and the more magnetic the re¬ 
sulting ground-state whereas the more positive T, the 
more antisymmetric the spin representation and the more 
Fermi-liquid like the ground-state. To avoid biasing the 
physics, we consider a grand-canonical ensemble of rep¬ 
resentations defined by the partition function with indef¬ 
inite asymmetry F, 

Z = Tr[PQ„e-^"] = Y,Tr[PQ,,Yoe-^% (12) 
Yo 

where now we can identify PQo with the large-A/" general¬ 
ization of the Gutzwiller projection operator introduced 
in Eq. 

/ rif) 

ll.^eiOi{nBi+nFi-Qo) ^ ( 13 ) 

27r 

This procedure will enable the ensemble to explore 
the lowest energy configurations. Another motivation to 
work with the constraint that fixes only the total number 
of boxes of the representation is the fact that the asym¬ 
metry of the representation appears only in the large-A" 
limit, so by letting Y run free provides an unbiased way 
to take the limit A ^ 2, as schematically shown in Fig.|^ 

C. New features of this work 

The large-A limit we now develop places all L-shaped 
representations with a given number of boxes on the same 
footing, and the asymmetry of the representation can be 
thought of as a variational parameter. The character 
of the representation (bosonic, fermionic or mixed) will 
now be decided by the energetics of the problem. This 
will permit us to explore the phase diagram of systems 
as heavy fermions, in which the character of the spin 
changes from fermionic in the HFL phase towards bosonic 
in the AFM region. 

Fig. [^illustrates schematically the evolution of the en¬ 
ergy landscape (energy as a function of gi?, the number of 
fermions in the representation), for three different values 
of Tk/Trkky- For small Tk/Trkky the energy land¬ 
scape has a minima for qp = 0, which means that the 
system prefers to have a bosonic spin representation and 
possibly develops magnetic order. Analogously, for large 
Tk ITrkky ^ the energy landscape has a minimum for 



Tk/Trkky 

FIG. 4. Schematic representation of the evolution of the en¬ 
ergy landscape as a function of for different values of the 
ratio TkITrkky ■ 


the maximum value of g^, indicating a purely fermionic 
representation, which would possibly lead to the devel¬ 
opment of a heavy Fermi liquid. For intermediate values 
of Tr/Trkky we find that the representation is mixed, 
with an energy minima developing at an intermediate 
value of qf so that the minima is a saddle point as a 
function of gi?. We are thus able to identify two classes 
of solution: 

- Type I minima, in which the free energy is mini¬ 
mized by a purely bosonic or a purely fermionic 
representation, indicating that the original super- 
symmetry of the spin is severely broken. In this 
case the results of a purely bosonic or fermionic 
representations are recovered; 

- Type II minima, in which mixed representations 
are energetically favorable. These kinds of min¬ 
ima are candidate representations for a two-fluid 
picture of heavy fermions. Since the fermionic and 
bosonic components of the spin fluid acquire the 
same chemical potential, this opens up the possibil¬ 
ity of a new kind of zero mode: a Goldstino, arising 
from the zero energy cost of rotations between the 
fermionic and bosonic spin fluid. 

Other new features of this work are: 

- Symplectic Spins: we generalize the rotation group 
of the spin from SU{2) to SP{N) for a large-A 
treatmenfP^^^. This generalization guarantees that 
all the components of the generalized spin properly 
invert under time reversal which confers various ad¬ 
vantages. In particular, it allows the description of 
geometrically frustrated magnetisnJ^ and it per¬ 
mits the exploration of singlet superconductivity 
within the large-A framework. 

- Spatially inhomogeneous representations: We ex¬ 
plore the possibility of solutions which sponta¬ 
neously develop Kondo, or magnetic character at 
different sites. Here we are motivated by the 
partially ordered phase verified experimentally in 
CePdAl. Fig. represents this kind of solution 
schematically in a frustrated triangular geometry: 
one of the local moments in the triangle develops 
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fermionic character, forming a singlet with elec¬ 
trons in the conduction sea, while the other two 
local moments have a bosonic representation, form¬ 
ing an antiferromagnetic bond. 



Local moment 
Conduction electron 


spin 


FIG. 5. Schematic representation of the inhomogeneous solu¬ 
tion for the frustrated triangular geometry. 


This paper is organized as follows: in Section [H] we de¬ 
fine the supersymmetric-symplectic spin, identifying the 
gauge group under which it is invariant and discuss the 
Casimir for a given irreducible representation. Details on 
the derivations are given in Appendix In Section O 


we introduce the general path integral formalism which 
leads to a mean field free energy in the large-A" limit. 
We apply this formalism to a two-impurity model in Sec¬ 
tion [IV| and details on the calculations are given in Ap¬ 
pendices 00 and In Section |V| we apply the same 
formalism to a three-impurity model. We conclude and 


discuss the open questions in Section VI 


II. SUPERSYMMETRIC-SYMPLECTIC SPINS: 
DEFINITIONS AND PROPERTIES 

We start defining the supersymmetric-symplectic spin: 

SaH = fifp - + bihn - aPb'_pb-a, (14) 

where /j^, and are respectively, fermionic and 

bosonic creation and annihilation operators, with indexes 
a, P = {±1, ±2,..., ±A/2}, and a = sgn{a), P = sgn{p). 
In this form the invertion of the spin under time reversal 
is made explicit. 

We can write the spin operator more concisely as: 


Saf} = «'l70«'/3 = (15) 

by introducing the four component spinor Tq, = 




a 


/ /« \ 

^f-a 


(16) 


\ / 

which carries the explicit spin index a, and has an im¬ 
plicit super-index which runs from 1 to 4 related to the 


supersymmetric and particle-hole character of its entries; 
and the matrix 


/ 1 0 0 0 \ 
0 10 0 
0 0 10 
Vo 0 0 -1/ 


(17) 


We shall follow the convention that super-indices are sup¬ 
pressed and fully contracted with one-another in our for¬ 
mulae, unless otherwise stated. 

The supersymmetric-symplectic spin defined in Eq. \TE\ 
commutes with the following operator bilinears and the 
respective conjugates: 

flF = 

a 

flB = '^b^ba, 

a 

'tp = ^faf-a, (18) 

a>0 

OL 

fi = y^^gfab-q. 


These are therefore generators of the symmetry group of 
the supersymmetric-symplectic spin. Note that t/i and fj 
are present in the SP{N) but not in the SU(N) gener¬ 
alization of the supersymmetric spiiP^. 

We can rewrite these operators in the form of Hubbard 
operatorJ^ as follows: 


^00 

X__ 

^+- 

^+0 

^-0 


{flF - Ub) 

2 

{N — ftp — bb) 


y, x_+ = V), 

71 ’ ^"+"71 

fj _ 

^ 2 ’ ^ 2 ’ 


(19) 

( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 


where Aqo, X±± and V±=p are bosonic Hubbard opera¬ 
tors, while X±o and Xo± are fermionic Hubbard opera¬ 
tors. In this form the algebra that these operators follow 
can be concisely written as: 


^cd]^ — ^ad^bc A ^cb^ad-) 

where the anticommutator, (+), is used when both Hub¬ 
bard operators are fermionic, while the commutator, (—), 
otherwise. The Hubbard algebrgP^ above defines the 
SU{2\1) supergroup!^. One can explicitly see the SU{2) 
subgroup generated by the isospin operators 


Ti = (X+_+X_+)/2 , (26) 

T2 = -i(X+_-X_+)/2 , (27) 

T3 = (X++-X__)/2 , (28) 
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which follow the commutation relation: 

(29) 

and are related to the rotation of the fermionic compo¬ 
nents of while Xqq defines the generator for the U{1) 
subgroup associated with the bosons. 

Note that one can perform a super-rotation g taking 
^ which leaves the spin invariant: 

^'l7o^'/3 ^ (30) 

if g satisfies 

g^iog = 70- (31) 

The most general transformation g can be constructed 
by exponentiation of the generators of the group listed 
above. An explicit form for g is given in Appendix ]A 2| 


A. The Casimir 


This form is sim ilar to that found for the SU(N) case (see 
Appendix |A4| for details of the derivation). Note that in 
the lowest weight state |To), where V^ITq) = = 

^|To) =0, (corresponding to no pairs and to minimizing 
the number of fermions) = 2{nB ^nF){N — np 

We can also write in terms of the SU{2\1) Casimir 
of the Hubbard operators via the identity: 

{N^ - S2)/4 = XapXpa - [XaO, ^Oa] - (Xoo)t (37) 


with an implied summation over the repeated indices 
(a,/3 = ± (See Appendix A5). 


III. THE FORMALISM 

The main topic of interest in our work is the class of 
Kondo-Heisenberg models involving SP{N) spins, inter¬ 
acting via an additional nearest neighbor antiferromag¬ 
netic Heisenberg exchange. These models are written as 


To uniquely characterize an irreducible representation 
of SP{N), in principle one needs to define r = N/2 — 1 
Casimirs, where r is the rank of the group (the dimension 
of the Cartan sub-algebra). In terms of Young tableaux, 
one can understand these parameters as the number of 
boxes in each row of the tableau (the maximum number 
of rows in the Young tableau for SP{N) is N/2). In this 
work we restrict our attention to L-shaped tableaux, so 
we have the extra information that all the rows below 
the first one have no more than a single box. This re¬ 
duces the number of parameters required to define the 
representation to two: Q, the total number of boxes, and 
Y, the asymmetry of the representation, as discussed in 
Coleman et al.^^. 

From the second Casimir we can identify the quantities 
Q and Y. From NwachukrP^, we can deduce that for an 
L-shaped representation in SP{N) the second Casimir 
can be written in terms of the width w of the first row 
and the height h of the column in the tableau as: 

C 2 = 2(rc + h){N Yw-h)Y4.{h- N/2) - 2, (32) 
and identifying Q = w Y h — 1 and Y = h — w, we have 


H = F + + E (38) 

3 ihj) 


where 


Hc = Y^ ek4„Cka (39) 

ko! 

describes a conduction band of electrons of dispersion Ck, 
where creates a conduction electron of momentum k, 
spin index a G {±1, ±2,..., YN/2}. The term 

HK{j) = P^Sa^{j)sp^{j) (40) 

a/3 


describes the Kondo interaction at site j, where Sag{j) = 
^ia^i/3 defines the local moment operator as in Eq. 


15 


and the conduction electron spin operators are also writ¬ 
ten in symplectic form: 


SaffU) = c]^Cjis - aj3c]_Fj-c<- ( 41 ) 

The final term describes the Heisenberg interaction be¬ 
tween spins at sites i and j, given by 


C2 = 2Q{NyI-Y), (33) 

the details of this derivation are shown in Appendix ] A 3| 
In terms of operators, the second Casimir can be writ¬ 
ten as the magnitude of the spin: 

= (34) 

a/3 

= 2Q{N + 1-Y) 


where; 


Q = hp + n-B, 

Y = hp - hs + 1 + 


Yip - + 27?t?7 

0 ■ 


(35) 

(36) 


HniiJ) = ^'^Sap{i)Spa{j)- (42) 

a/3 

To display the supersymmetric gauge character of the 
interactions we now rearrange the order of the opera¬ 
tors. First, using the properties of the symplectic spin, 
^PSaf 3 = —S-g-a^ the two parts of the electron spin 
operator in the Kondo interaction are folded into one as 
follows: 

Hk^J) = Sai3{j){c^jpCja — dLp(^-_^Cj-g) 

a/3 

= ^ E ^«(j)^'/3(i)c]^9«- (43) 

a/3 






7 


Next, we super-commute the ^ field to the right-hand 
side of the interaction, rewriting the interaction in terms 
of a supertrace. 


form like and under super-rotations. 

In a similar fashion, we rewrite the Heisenberg inter¬ 
action as 


HkU) = 


N 


al3 


(44) 


where we define the supertrace as Str[H] = An A 22 — 
Ass — An and we have used the property that the dot 
product of two super-spinors and y can be rewritten as 
a supertrace of their outer-product [4>x] = —Str[x<h]. No¬ 
tice that Uj = and Uj = (^Cjp^jp) are four com¬ 

ponent column and row spinors that respectively trans¬ 


= “W • (45) 

Notice that the object Uij = {^ijs^jjs) is an outer- 
product of two super-spinors, forming a four-by-four ten¬ 
sor in superspace that transforms as Uij QiUijgJ^ un¬ 
der super-rotations. 


With these manipulations, the Kondo-Heisenberg model can be written as: 


H = Hr- 


2Jk 

N 


^Str 

j,al3 


jg) 


Jh 

N 


^ Str . 

ii,j)a,/3 


(46) 


Notice that the invariance property of the supertrace Str[^T^“^] = Str[T] guarantees that these interactions are 
gauge-invariant under the local super-rotations The factorized forms of the interactions are convenient 

for Hubbard-Str at onovich transformations. 


The constraint fixing the total number of bosons plus 
fermions at each site npi^riBi = Qo can also be written 
in terms of the spinors 


where 


npi 


■nBi = 7: 


/ 1 0 0 0 \ 

0-100 

0 0 10 

\ 0 0 0 -1 / 


(47) 


(48) 


is the measure of integration over the canonical c, / and 
b fields and 


'DX = '[[dXj (52) 

3 

is the measure of integration over the constraint. 

The components of the action are: 

rg 

Sc= + Ck)Cka, (53) 

•^0 ka 

the conduction electron part of the action; 


We can now write the partition function as a functional 
integral over the constraint field A and the spin carrying 
boson and fermion fields 

describes the Berry phase, with the constraint ubj + 
VXDiie~^ (49) ~ imposed via the introduction of the Lagrange 

’ multiplier Xj at each site, while 



dr 


i ^ VQo (54) 




where 


S = Sc + Ss^Sk + Sh^ 


while 


T>ij. = D[c, f,b], 


(50) 


( 51 ) 



(55) 


are the Kondo and Heisenberg parts of the action. 


Inside the path integral, we can now carry out a Hubbard-Str at onovich transformation of the interactions. The 



Kondo part of the interaction is factorized as follows 

2Jk 


HkU) = Estr 

al3 

H'kU) = E Str r,- + y,. (c, 


= Ef(^i^. 


^ja 




N 

2Jk 


Tr[V,V,] 


(56) 


In the last step we have absorbed the minus sign associated with the anticommutation of ct^ and (Vj'^ja), likewise 
Cja and (^jaVj). These Hubbard-Stratonovich fields Vj and Vj = Vjjo are four-component spinors 





~ ('^i5 dj , <pj , ^j). 


(57) 


Here Vi and di are complex fields related to the hybridization between f-fermions and c-electrons and the development 
of superconductivity by the formation of pairs between f-fermions and c-electrons, respectively. The parameters 
and are complex Grassmann numbers, the first related to the hybridization between b-bosons and c-electrons and 
the second related to the development os pairs formed between b-bosons and c-fermions. 

In a similar fashion, the Heisenberg term decouples as: 


^ ^ ^ Str {^ja'^ia) T (^ia^ja) A^jj -|- j Str I^A^jfA^jf] 

OL 

= + ^jaAij'^ia] + Tstr [AijAij] (58) 

a 

where Aij is a four-by-four matrix and its conjugate is defined as Aij = yoA^yo- The structure of the matrix is as 
follows 


where the diagonal block-matrices composed by c- 
numbers whereas the off-diagonal block matrices As and 
A^ are composed by Grassmanians. The internal struc¬ 
ture of these blocks is given by 

= 5(). (») 

= ( X %, ) ■ <®> 

where the remaining matrix is defined through the rela¬ 
tion {As)ij = (T 2 ,{As)ji. The matrix Aij can be thought 
of as supersymmetric RVB field. The components fields 
tij and pij promote hopping and pairing amongst the 
f-fermions in different sites and the complex fields qij 
and gij promote hopping and magnetic bond formation 
between the b-bosons. The Grassmannian parameters 
^ij are hopping amplitudes that transmute bosons into 


fermions and vice versa, while the Grassmannian ampli¬ 
tudes Pij describe pairing between bosons and fermions 
at different sites. 


The partition function now reads: 


z = j V[\V,A]ViJLe-^\ 

(63) 

where 


S' = + S'h, 

(64) 

where the primes denote the actions of the 
Stratonovich factorized interactions 

Hubbard- 

S'k = [ dr 
do 


3 

s'j= [ dTYH'nihj), 

{ij) 

(65) 


and the integral over P[H, A] indicates the integral over 
all the fluctuating fields introduced by the Hubbard- 
Stratonovich transformations. Here c-number fields are 
represented by the latin letters (i;, (i,p, t, g, ^), while the 
Grassmannian fields are represented by the Greek letters 
(0, (^,y,/i). Grassmannian fields are introduced in order 
to decouple terms with fermionic bilinears. Note that the 
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Kondo and Heisenberg parts of the action are invariant 
under the transformation ^ dj^jcr if the fluctuating 
field matrices transform accordingly as 


V 9 j Vj 

(66) 

^ij ^ 

(67) 


Now we move to the discussion of the implementation 
of the constraint by fixing Q = Qq. The constraint can be 
imposed as a projection operator in each site j, written 


as a delta function: 

Pq. = ^jPL = - Co). (68) 

In order to treat the bosonic and fermionic components 
of the spin in the grand canonical ensemble, we split the 
constraint into two terms as follows: 

Qq 

Pqo ^ ~ QFj)S{nBj -Qo + Qpj)- (69) 

QFj=0 

The constraint fixes the total number of bosons and 
fermions to Qo, leaving the asymmetry of the represen¬ 
tation free to adjust according to the energetics of the 
problem. 


This constraint can be implemented in the path integral as a Dirac delta function in its integral form: 

Pqo = 1 

T>[Qf, A] = HE dXpjdXBji 

3 Qfj 


where 



XpjiriFj — Qfj) + XBjiriBj — Qo + Qfj) , 


(71) 


and the constraint fields, Xpj and Xbj are integrated along the imaginary axis. From these considerations Ss can be 
rewritten as: 


Ss ^ S's = f X] f E - ^FjQpj - ^BjiQo - Qfj)^ , 


where now 


j \ a 




^ Xpj 0 0 0 

0 ~XFj 0 0 

0 0 Xbj 0 

Vo 0 0 ~XBj 


and the partition function is now written as: 


Z = 




where 


S' = S, + S's + S'f + S'h. 


(72) 


(73) 


(74) 


(75) 


The new feature of this action, is the appearance of the 
term Qf, which tunes the bosonic/fermionic character of 
the representation. In the large N limit, we will be able 
to replace the discrete measure of integration over Qf by 
a continuous measure 



{N ^oo). 


(76) 


In the large N limit, we anticipate that the functional 


integral is given by the saddle point value of the effective 
action, so the large N approximation is then given by the 
exponential of the effective action 




(77) 


where 


e 



e 


-s' 


(78) 
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is the integral over the canonical spin-carrying fermions 
and bosons in the presence of fixed Qp A^, Ai?, V and A 
with the conditions that Sef / be stationary with respect 
to each of its fields. We may anticipate two classes of 
mean-field solution 


1. Type I solutions, in which the minimum of the ef¬ 
fective action occurs at the extremum of the sum¬ 
mation over Qpji he Qpj = Qo oi* Qfj = 0 , corre¬ 
sponding to the fully fermionic or bosonic solutions. 


2 . Type II solutions, in which the minimum of the 
effective action occurs at some intermediate value 
of 0 < Qpj < Qo- Since the action is stationary 
with respect to variations in Qpji this implies that 


SQfj 


= 0 = Xpj — Xpj 


(79) 


SO that in this phase, the chemical potential of the 
bosonic and fermionic spinons are equal. 


\Bj=XFj. (80) 

This equality of chemical potentials allows us to 
consider these solutions as two fluid solutions. 


connections to the physics of heavy fermion systems. The 
model consists of two local moments interacting among 
themselves by a Heisenberg coupling Jp and interact¬ 
ing with its respective bath of conduction electrons by a 
Kondo coupling Jk- 



^ Lead 2 j 


Jk 


FIG. 6. Schematic representation of the two-impurity model. 


The Hamiltonian is written as: 

H = Hc + JkYI (0)^a/3a + SiapS2/3a{83) 

a,a/3 af3 


At the saddle points, we can set all fermionic components 
of V and A to zero. These terms only contribute to the 
fluctuations about the mean field theory, so we have 


and 


v;- ^ w 



A^- 


Ap 0 \ 

/ ij 


(81) 


(82) 


Note that the fermionic and bosonic parts of the action 
decouple since all the matrices in the action now have 
the blocks linking the fermionic and bosonic subspaces 
equal to zero. Now it is possible to solve the fermionic 
and bosonic problems separately, imposing the constraint 
Qfj^QB j = Qo to the solution in the end of the calcula¬ 
tion. Note that this provides a picture of two asymptot¬ 
ically independent fluids, bosonic and fermionic, in the 
large-A limit and that the introduction of fluctuations 
will provide interactions between them. 

Now, as a first exploration of this idea we illustrate the 
formalism with two simple examples: a two-impurity and 
a frustrated three-impurity model and show that there 
are stable mean field solutions with mixed representa¬ 
tions, as well as with purely bosonic and purely fermionic 
representations. 


IV. THE TWO IMPURITY MODEL 


where He = is the conduction electron 

Hamiltonian, a = {1,2} is the lead and local moment 
index, k the momentum and a the spin index, which as¬ 
sume values cr = {±l,±2,..., ±N/2} in the large-N limit. 
Here 5 ^( 0 ) is the spin density of conduction electrons 
at the site which is connected to the local moment spin 
Sa- Introducing the supersymmetric-symplectic spin, the 
Hamiltonian can be written in the large-V limit in terms 
of fermionic and bosonic operators as: 

F = F, - ^ ^ Str {c,0^a0)] 

a,a/3 

- ^ E Str [(^'l/3^2;3) (^' 2 a^la)] • (84) 

a/3 


where Ter = Tj. 70 , as defined in Eq. 

We can now apply the formalism introduced in the 
previous section. We perform a Hubbard-Stratonovich 
transformation to decouple the interacting terms in the 
Hamiltonian by introducing fluctuating fields. Within 
a static mean field solution the fermionic and bosonic 
problems decouple and are effectively linked only by the 
constraint; we can now factor the partition function as: 

Z[qp] = Zp{qp)ZB{qo - qF), (85) 

where we have defined go = qF = Qf/N- We 

see solutions where the Free energy is minimized with 
respect to qp. 

The fermionic part of the partition function reads: 


16 


As a first application of the supersymmetric-symplectic 
spin, we study a minimal model that allows one to make 





o~^F 


V^lF=V[cJ] 
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Sf= S, + [^ dr [flidr + \F)faa 
a,cr 


+E a^akcr H“ h.C.^ 

k 

+/3Af y ^ - 2pN\FqF, (86) 


where we already dropped the terms in Pa and ta^ since it 
can be shown that these do not contribute to the saddle 
point solution. Also, the f-operators can be redefined 
to eliminate the pairing term between c- and f-operators 
from the Hamiltonian: 


faa ^ fa 


Vafaa + daafl_^ 

VkP + KP 


(87) 


Under these considerations, the fermionic part of the 
solution reduces to two decoupled impurity problems. 
Taking Va = v to be site independent, integrating out 
the conduction electrons in each lead and transforming 
from imaginary time to Matsubara frequencies (see de¬ 
tails in Appendix [B|) : 


Sb 


I 




■l5l^ 


+/3N— - 2/3NXB{qB + 1/2), 

Jh 


where 


Lb = 


dr + Xb g 
9 —dr + Xb 


(95) 


(96) 


= 


f ha 

V ^^l-a 


(97) 


Here we already dropped the fluctuating field q since the 
saddle point solution results in g = 0. Integrating out 
the bosons and summing over Matsubara frequencies (see 
Appendix]^, in the zero temperature limit, the free en¬ 
ergy is given by: 

^ - \g? + ^ - 2AB(te + 1/2). (98) 

Minimizing the free energy with respect to g and \b 
one finds: 


^ ^ T T n)faerifyn) 

naa 

N\v\‘^ 

+2/3—-- 2f3NXFqF^ ( 88 ) 

Jk 

Tn =Te{D-\iujn\)sgn{iujn)^ T = Fpo\v\‘^, (89) 


Ab = Jh{qb + 1/2), = JjjqBiQB + 1), (99) 

so the the bosonic free energy can be written as: 

= -JniqB + 1 / 2 )^, ( 100 ) 


where po is a constant DOS, D is the bandwidth and 0(x) ^ constant term, 

is a Heaviside step function. Summing over Matsubara 
frequencies, in the limit T ^ 0, the free energy has the 

r A. Analysis of the free energy 


N 



CfIii 


_Y 

J\ ’ 


where we define 


(90) 


= Ab + 'i^r, (91) 

and the Kondo temperature 

Tk = (92) 

In the large-A limit the partition function is dominated 
by the saddle point. Minimizing the free energy with 
respect to one finds 


= (93) 


and substituting back into the fermionic free energy: 

^ =--T/f sin(7rgF). (94) 

iV TT 

The bosonic part of the partition function can be con¬ 
cisely written as: 

Zb = J VpBe ^~^^, T^Pb = 9, Ab] 


Making explicit use of the constraint condition, gB + 
Qb = Qo^ the total free energy can be written as: 

-f— = --AsihnqF) - {qo - qF + l/2)^ (101) 

Jb/V TT 

where A = ^ and the free energy is given in units 
of Jh- For each value of A and go the representation 
was determined by the minimization of the free energy 
with respect to gB and the result is plotted in the rep¬ 
resentation diagram of Fig. 0 In case the free energy 
is minimized for gB = go the phase is purely fermionic, 
meaning that a completely antisymmetric representation 
is favored. Analogously, for gB = 0 (or gB = go) the 
phase is purely bosonic, and a symmetric representation 
is more appropriate. Solutions with 0 < gB < go are solu¬ 
tions in which both bosons and fermions coexist, which 
we call a mixed phase and label as (F + B) in Fig. 
Note that for a fixed value of go, as the ratio TkIJh 
is increased the spins tend to develop fermionic charac¬ 
ter. Also, for a fixed value of TkIJhi increasing 1/go 
(or reducing go, which is equivalent to decreasing the 
magnitude of the spin) the spin representation also tends 
towards a fermionic representation. 
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nF=qo 

nB =0 


TkUh 


FIG. 7. Representation diagram for the two impurity model 
indicating the most favorable representations as a function of 
A = ^ and l/qo- A purely fermionic representation phase 
(red) is labeled hy F, a pure bosonic representation phase 
(blue) is labeled as B and the mixed representation phase 
(intermediate colors) is labeled by F B. 


The dashed line between F and F-\-B regions indicates 
a second order phase transition and can be determined 
from the condition: 


which leads to: 


OF 

dqF 


= 0 , 

qF=qo 


A = 


1 

2 cos(7rgo) * 


( 102 ) 


(103) 


The continuous line represents a first order phase tran¬ 
sition. The line between the purely fermionic and purely 
bosonic representations is determined by: 


Fjqr go] _ F[qF 0] 
N ~ N 


(104) 


which gives the condition: 


^ ^ 7rgo(go + 1) 

2sin(7rgo) 


(105) 


The first order line between the phases B and F B 
cannot be computed analytically and was determined nu¬ 
merically. 

Throughout the mixed phase we have Xf = A^, so 
both fluids have the same chemical potential. This is 
related to the presence of the Type II minima of the free 
energy (see discussion in the introduction), with a saddle 
point that allows the coexistence of bosons and fermions 
and the interchange of one into another at no energy cost. 
In particular, at the second order “phase transition” line 


discussed above, when qF = Qo — S and qs = ^ with 
^ 0, one can check explicitly from Eg. [M| and that 
the condition Xf = Xb gives the same condition that 
defines the second order line in Eg. |1Q3[ 


B. Fluctuations of the local fermionic fields 


We now analyze the effects of fluctuations of the lo¬ 
cal fermionic fields. In the previous section we in¬ 
troduced the time dependent fields and which 
allow us to decouple the terms 

and (d&I_a4^^)(4cak/3&a-/3) in the action, respectively. 
These fields do not acquire an expectation value, but flue- 
tuate around zero. The partition function can now be 
written in terms of the saddle point solution determined 
in the former subsection times and the new con¬ 
tributions to the partition function due to the presence 
of the fluctuating fields Sef) and that we take to be site 
independent. Eocusing first in the Sc/) field: 



fo dr |^Eaka(^Lcaka(5(/)+/l.C.)+^Ai^ 


(.106) 


Expanding to second order in be/) we can identify the 
propagator for the fluctuating field Sc/): 


= 2N 


Xeb 


1 

Jk 


(107) 


with 


XebiX^r) 


, ' ' . baa(iwm+iain) 



where Gb{iUn) is the bosonic propagator and Gkk'(^<^n) 
is the full c-electron propagator. 

We now evaluate XebJ^) (details of the calculation can 
be found in Appendix |M. One interesting region for the 
analysis of D(j){iuJr) is the second order transition line, 
where the energy levels of the bosons and fermions are 
egual. In the infinite bandwidth limit we find: 


Xcb{(^ - iS) - ^ = poujRe 
Jk 


log(A -b iA) 
cj + iA 


Polo‘s 

A 2 


log(A — uj F ib). (108) 


Note that at zero freguency, Xc6(0) ~ ^ = 0, so 
[D(f){0)]~^ = 0, and the propagator for the fermionic hy¬ 
bridization field 0 diverges at zero freguency at the sec¬ 
ond order phase transition, indicating the presence of a 
fermionic zero mode. Also, there is a gap of magnitude 
egual to X = Xf = Xb with a continuum that goes up 
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to the bandwidth. This gap is always present in the 2- 
impurity model since \ = = Jh aways finite at 

the transition. For a Kondo-Heisenberg model in the lat¬ 
tice, the bosonic level will acquire a dispersion and when 
magnetic order sets in it will be gapless at some points 
in the Brillouin zone. In that case the spectrum for the 
fermionic hybridization field is expected to have a con¬ 
tinuum of excitations, which can potentially lead to non 
Fermi liquid behavior. 



(ji) 


V. FRUSTRATION IN THE THREE IMPURITY 
MODEL 

As a second application of the supersymmetric- 
symplectic spin, we study a minimal model that brings 
in the issue of geometric frustration into play. The model 
consists of three local moments interacting among them¬ 
selves by an antiferromagnetic Heisenberg coupling Jh 
and interacting with its respective bath of conduction 
electrons by a Kondo coupling Jk^ as depicted in Fig.[^ 
We are motivated to look at this problem by experiments 
in CePdAl in which the equivalent Ce sites spontaneously 
develop a state in which one third are paramagnetic and 
the other two thirds are magnetically ordered^, explor¬ 
ing the ability of the symplectic representation of the spin 
to describe frustrated systems. 

The Hamiltonian is written as: 

H = H, + JkY. 

^aa p{0)SaPa Jh ^ ^ ^aQ:/3^0+1/30:; 

a,a/3 a,a/3 

( 111 ) 


FIG. 8. Plot of the imaginary time of D{uj)^ the 6(j) propa¬ 
gator. The parameters used in this plot were a solution of 
the mean field theory at a specific point of the second or¬ 
der phase transition: A = 1.57 and 1/qo = 2.5, which gives 
Af = As = A = 1.01 and A = 2.97 (in units of Jh)- 


For the second fermionic mode a similar calculation 
follows, where now: 


= 2N 


\cb(J'^r) j 

Jk 


(109) 


with 


, - ' . bao(icjOm-iCjOn) 



and we can write, at the second order transition line: 
Xc6(w -i 5 ) - A = 2\)Re 

■JK 

Here we note that there is no zero mode for this 
fermionic field at the transition, and a continuum starts 
at a finite X = This is related to the choice 

we made to implement the constraint, which is not in¬ 
variant under all transformations that leave the spin in¬ 
variant. In particular, it is not invariant under transfor¬ 
mations of the form qb (see Appendix [A|), generated by 
the operators 77 ^ and 77 . 


log(-A -I- iA) 
uj — 2 A “h iA 


where now a = {1,2,3} is the lead and local moment 
index with periodic boundary conditions. As in the pre¬ 
vious section. He is the conduction electron Hamiltonian, 
5a ( 0 ) is the spin density of conduction electrons at the 
site that is connected to the local moment spin Sa- 
Introducing the supersymmetric-symplectic spin from 
Eq. into the Hamiltonian, this can be written in the 
large-A limit as: 

F = Fc - ^ E Str [(^aacL) (Ca/jta/j)] 

a,a/3 

- ^ E Str [{'ifaP^a+10) (^-a+la^aa)] • (112) 

a,a/3 

with Tacr defined in Eq. [T6| 



FIG. 9. Schematic representation of the frustrated three- 
impurity model. 
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Proceeding as in the previous section, within a path in¬ 
tegral formalism we introduce fluctuating fields in oder to 
decouple the quartic terms in the action and impose the 
constraint by the introduction of a delta function in the 
integral form. Within a static saddle point solution the 
problem decouples in to a bosonic and a fermionic part, 
effectively linked by the constraint. In this section we 
are going to leave the representation of the spin and the 
mean field parameters to be determined independently in 
each site. Omitting the details (similar to the previous 
section), the partition function in the large N limit can 
be written as 


Z = ZfZb, (113) 

with the understanding that the partition function is to 
be stationary with respect to the qpj at the three sites. 

The fermionic part of the partition function can be 
written as: 

Zf = J ViiFe~^^, VfiF = ^[c, /] 

r/3 

Sf= Sc + rfr V [faaidr + >^Fa)faa 

Jo 

+ a^akcT T h.C.^ 

k 

+/37V^t^-/37V^AFa?Fa, (114) 

a a 

where as in the previous section, we assume that the only 
fluctuating field that acquires a finite value at the saddle 
point solution is Va- In this case the fermionic part of the 
solution reduces to three decoupled impurity problems, 
with the same solution as the previous section, now for 3 
leads: 

^ = sin(7rgFa)- (115) 

a 

The bosonic part of the partition function reads: 

= Jv[b]e-^^, (116) 

where 

Sb = 

a 

- PN ^ ^BaiQBa + 1/2), (H^) 

a 

where we denote QBa = Qo — QFa and 

( tubi A^i A^3 \ 

^^Bl Ab 2 I 5 (118) 

Ab3 A^2 / 


with 


BflBa = 


dr + ^Ba 

0 


0 

— dr + A^a ) ’ 


(119) 


A^a = 


and 


^ Ba = 


Qa 9a 

9a Qa 


\ 

^2(7 


( 120 ) 


^3cr 

V ^ 4-0 


( 121 ) 


/ 


Note that the trace term in the action now appears with 
a minus sign since it is related to the bosonic part of the 
super-trace introduced in Eq. 58 Here we define Jq = <73 


as the bosonic part of the original matrix 70 . 

The solution of the bosonic part of the partition func¬ 
tion is more involved if we allow the representations of 
the spin to be different in each lead. As a first solu¬ 
tion we consider the same representation on every site, 
and look for an homogeneous solution^ taking A^a ^ A^, 
Qa ^ Q and 9a ^ 9- Integrating out the bosons and 
summing over Matsubara frequencies, in the zero tem¬ 
perature limit, the free energy can be written as: 

Fb 


= (Ab + 2g) + - qY 3g‘^ 

H —~ ^^b{qb + 1 / 2 )- ( 122 ) 

Jh 

Minimizing the free energy with respect to g, g and \b 
one finds the saddle point free energy: 


N 


3Jh 


{qb + 1/2)^ 


(123) 


up to a constant term. The total free energy for the 
homogeneous solution can be written, already making 
explicit use of the constraint condition qf F Qb = Qo^ as: 

= - lAsin(7rgF) - ^{qo - Qf + 1/2)^ (124) 
JhI\ tt Z 

where again A = ^ and the free energy is given in units 
of Jh- Note that this is functionally the same as the 
2-impurity model up to an overall factor of 3/2, and as 
a consequence the representation diagram determining 
the most favorable representation for the spin within an 
homogeneous solution will be identical to the 2 -impurity 
case. 

Now we move on to investigate solutions which spon¬ 
taneously develop different representations in each site, 
which we refer to as inhomogeneous representations. Due 
to frustration, we expect that it is energetically favorable 
for one of the spins to be in a fermionic representation, 
essentially disconnected from the other two spins with 
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Inhomogeneous 


nB=qo 



np-qo 

nB=0 


TkUh 

FIG. 10. Representation diagram for the three impurity 
model indicating the most favorable representations as a func¬ 
tion oi A — ^ and 1/go- The color code is the same as in 
Fig. The hashed area represents the region of the diagram 
in which an inhomogeneous solution is energetically favorable. 


a bosonic or mixed representation, forming an antiferro¬ 
magnetic bond. We assume one of the spins to always 
have a fermionic representation and let the representa¬ 
tion of the two other spins to be selected as the one that 
minimizes the total energy. 

Now the problem reduces to a single impurity prob¬ 
lem with a purely fermionic representation plus the two- 
impurity problem solved in the previous section. The free 
energy for the inhomogeneous solution reads: 

1^(7 =-l^sin(7rgo) --^sin(7rgF) (125) 

JhAI tt tt 

- {qo -qF + 1/2)^, 

in units of Jh, where A = ^. Again, the representation 
diagram will be the same as before, but now we com¬ 
pare the free energies of the homogeneous and inhomoge¬ 
neous solutions for the 3-impurity problem. The hashed 
area in Fig. is the region of the diagram in which 
the inhomogeneous solution is more favorable. This re¬ 
sult provides a model for the situation which appears to 
occur in CePdAPl one third of the local moments in 
the frustrated Kagome lattice (formed by an assembly 
of corner-shared triangles) relieve the frustration by as¬ 
suming an antisymmetric character and forming a Kondo 
singlet with a conduction electron, while the other two 
thirds of the local moments assume a bosonic character, 
developing magnetic order, allowing a partially ordered 
phase to be formed, as depicted in Fig. 


In this work we introduced a new supersymmetric- 
spin representation for large-V treatments, based on 
a symplectic generalization of the spin operator. We 
have analyzed the properties of the supersymmetric- 
symplectic spin and its symmetries, identifying the super¬ 
group SU{2\1) of transformations under which the spin 
is invariant. 

We have proposed a new framework in the large-V 
limit, which allows the problem to sample different rep¬ 
resentations, selecting the one which lowers the energy in 
a given point in parameter space. This opens up the pos¬ 
sibility of describing the phase diagram of heavy fermions 
within a single approach that can capture the evolving 
character of the spin, at the same time that it offers a 
potential framework for the phenomenological two-fluid 
picture for heavy fermions. 

Applying this approach to two toy models, the two- 
impurity model and the frustrated three-impurity model, 
we have shown two new classes of mean-field theory that 
may be of interest in developing a unified description of 
heavy fermion systems. In particular, we find a mixed 
phase solution in which bosons and fermion coexist at 
each site, which points the way to a description of the 
coexistence of magnetism and Kondo effect. Also, we find 
stable inhomogenous solutions^ which may provide a basis 
for describing the partially ordered state in CePdAP^H^^. 

There are many open questions raised by our work. 
Firstly, within the mixed phase, we find evidence for a 
new kind of zero mode, a Goldstino, that results from the 
partial breaking of supersymmetry. This can be under¬ 
stood as a consequence of the fact that the Hamiltonian 
is invariant under super-rotations while at the same time, 
the mean-field solution breaks this rotational symmetry. 
Given a state |^) with a fermion assigned to the corner 
box, for example, one can rotate this state as follows: 

= (126) 

and find a new state which has same energy. These kinds 
of zero modes are only present when the chemical poten¬ 
tial of the bosons and fermions are equal, i.e in type II 
solutions. We note that in the presence of a Kondo effect, 
the f-electrons are charged, whereas b-bosons are neu¬ 
tral, so the fermionic Goldstino excitation is charged and 
may represent a zero energy valence fluctuation. Further 
work is needed to establish whether such zero modes are 
a truly physical excitation and whether they participate 
in anomalous inelastic scattering and the development of 
non-Fermi liquid behavior. 

Secondly, we would like to discuss the future appli¬ 
cation of this approach to the exploration of the phase 
diagram of the Kondo lattice. Extrapolating the results 
found in this work to the lattice we expect to find that 
for small ratios of TkIJh a bosonic representation is 
more favorable in which magnetic order will emerge once 
the bosons condense. Heavy Fermi liquid behavior is ex¬ 
pected to develop when there are fermions in the repre- 
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sentation, which can hybridize with the conduction elec¬ 
trons. If a mixed phase is stable in the lattice, heavy 
fermion behavior can coexist with magnetic order. In 
principle, this approach allows us to explore the different 
kinds of phase transitions as seen in YbRh 2 Siwithin 
a single approach: given that the fermions now can form a 
Fermi surface, the magnetic phase can emerge both from 
local moments within the bosonic representation or as 
an instability of the large Fermi surface. Another aspect 
of interest is the possibility of the description of super¬ 
conductivity; this phase is likely to develop in the mixed 
phase and its vicinity, as a consequence of a valence-bond 
kind of magnetism emerging from the fermionic antifer¬ 
romagnetic bonds. This is another interesting direction 
for future work. 
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Appendix A: Details on the properties of 
supersymmetric symplectic spin 


In this appendix we discuss a few properties of the 
supersymmetric symplectic spin in more detail. 


1. Number of spin components 


2. Symmetry transformations of the 
supersymmetric symplectic spin 


In the main text, a more concise form of the 
supersymmetric-symplectic spin is introduced in order to 
make the local symmetry manifestly clear: 

Sap = (A2) 


where. 


I fi \ 

y dib—fj j 


(A3) 


is a four-component spinor, and as defined in the main 
text 7 o = diag[l^ 1 , 1 , — 1 ], and a = sgn{a). 

Under a super-rotation of the spinors Tq, —the 
spin operator transforms as: 

Sa0 = 4'l7o«'/3 ^ '^ig^'yog'^0, (A4) 


SO for the spin to be invariant the transformation g should 
satisfy: 

9^09 = 70, (A5) 


which is essentially the unitarity condition to the trans¬ 
formation after taking appropriate care of the commuta¬ 
tivity of the bosons. 

The most general transformation g can be obtained 
by exponentiation of the generators of the algebra intro¬ 
duced in Eq. 18 Exponentiation of the even (or com¬ 
muting) part of the algebra gives: 


/ u 0 0 \ 

—fx 0 0 
0 0x0 ’ 

Vo 0 0 xj 


(A 6 ) 


When writing the components of the supersymmetric- 
symplectic spin as: 

Sap = flfp - a^'flpf-a + bibp - a^'blpb.a,{Al) 

note that we have in fact only N{N + l)/2 components 
(equal to the number of generators of the enlarged sym¬ 
metry group of the spin, here SP{N)). This can be 
checked by noticing that Sag = —a^S-g-a^ so the gen¬ 
erators are not all independent. Eor {3 > 0, we have 
Sag = —S-g-a, SO the generators with both indexes neg¬ 
ative are linearly dependent on the generators with both 
indexes positive. Eor a > 0, /3 < 0, a 7 ^ —[3 we have 
Sag = S-g-a, but uote that for the case of a = P we 
have the condition Sa,-a = Sa,-a^ which does not give a 
relation between the generators in the off diagonal. One 
can check that indeed Sa,-a is not linearly dependent 
on S-a,a' Within these considerations the number of 
independent generators is N{N + l)/2, as expected. 


where the parameters xx, v and x are complex num¬ 
bers satisfying |xxp + \v\‘^ = 1 and |xp = 1. Note the 
SU{2) and U{1) substructure of this transformation for 
the fermionic and bosonic parts of the spinor Ter, respec¬ 
tively. 

Exponentiating the odd (or anticommuting) part of the 
algebra we find: 


9A 

and 


9b 


/ 1 + ^ 0 
0 1 + ^ 
a 0 

V 0 —a 


—a 0 \ 

0 -a I 
1 - ^ 0 
0 1 - ^ / 


(A7) 



0 

0 


-/? 


0 Q -P \ 

1-7 ^ 0 

-p 1 + 7 0 _ 

0 0 1 + 7 / 


(A 8 ) 


where the parameters a and (3 are complex Grassmann 
numbers. 
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The most general transformation can be obtained by 
the composition of the three transformations above: 

9 =9EgA9B, (A9) 

which satisfies the condition g^^og = 70) as required. 


3. Derivation of the Casimir 

We now follow NwachukrP^ order to derive the 
second Casimir of SP{N), (where N is an even num¬ 
ber), given in Eq. Each irreducible representa¬ 

tion of SP{N) is characterized by the set of integers 
(/Ar/ 25 /Ar/ 2 -i 5 /i)- This set of numbers corresponds 
to the number of boxes in each row of the respective 
Young tableau, starting from the topmost and longest 
row, which has /tv /2 boxes. These numbers also cor¬ 
respond to the eigenvalues of the N/2 Cartan (diago¬ 
nal) generators in the highest state of the representation. 
EollowinJ^, if we define: 


Af — /i T N/2 T i, for i ^ 0 (AlO) 

\-i = —\i + Y, for i > 0 (All) 

Pi = N/2 + i, (A12) 

where we take /o = 0, the second Casimir is written: 

N/2 

C 2 = Y. (A13) 

i=-N/2 


Eor an L-shaped tableau with width w and height h, 
we have: 


fi = < 


w, i = f-, 

N 


1, f -l>i>f-/^ + l, 
0 , i < - h, 


(A14) 


therefore: 


\i — < 


i - N 

i 2 5 

i = -N 

i 2 5 


[wPN, 

-w, - 2 

1 + N/2 + ^ — 1^'^^w — h-hl, 


lo, 


\N 

I 2 

otherwise. 


-1 + N / 2-|-i, ~|"2-^1 — ^ — ~\^ - 


(A15) 


and 


\N, 

0 , 


Pi = \ 


N/2 + i, 
N/2 + i, 
.0, 


y — (A 16 ) 

-If-l|<i<-|f 

otherwise, 


so we can perform the sum: 


C 2 = + ANw + (A17) 

Jv/ 2-1 

+ Y {l + N + 2i) 

i=N/2-h-\-l 

-N/2+h-l 

+ Y 

i=-Ar/2+l 

= 2(^w -|- h)(Y w — h) T 4(h — Y/2) — 2, 


where the sums were evaluated as sums of arithmetic 
progressions. This is Eq.j^in the main text. Identifying 
Q = w P h — 1 and Y = h — re, we have 


C 2 = 2{Q P 1)(Y - Y) -f- 2(Y + Q + 1 - Y) - 2 
= 2{Q + 1)(Y - Y) - 2Y + 2(Y + Q) 

= 2Q(Y + 1-Y), (A18) 


which is Eq. similar to the form discussed in Coleman 
et al.l^ in the case of an SU(N) generalization of the 
supersymmetric spin. 


4. Operator form of the Casimir 


Now we relate the magnitude of the spin with the Casimir computed above. We start computing: 






I3a 5 


aP 


(A19) 
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with Soii 3 is defined in Eq. 14 


Taking the operator products and relabeling the summed indexes we find: 


+ blb/sblba + 2flfpblb„ 

- - a^bpblj.^ - 2apfy^blj.0), 


(A20) 


where we have used summation convention for the repeated spin indices. The first line of this expression is the Casimir 
for the SU{N) spin, while the second line introduces the additional cross-terms that result from the symplectic form 
of the spin. 

Expanding the first line, we have 

flfpflfa = nriN - np) + np 
babpb^ba =nB{N + ub) - ub 

fifpb^h^ = - np (A21) 

SO that 


flfpflfa + blbpb^i^ba + 2flfi3b^pba = {ub + nF){N + np - np) + 29"'9 - [np + np) (A22) 

Now expanding the additional cross terms, 

- aPfifisf-af-0 =nF- = np - 

-aPb'lbi3b^_^b-0 = ub - = np, 

-2apfif0blj.p = -2a'pflblj.0fp = -2v^v (A23) 


where we have used the fact that the bosonic pairs vanish pb-^bj^ = 0. Combining the additional cross terms, we 
have 


- aMlfpf-af-P - a^blbpbl^b-p - 2a^flfpbl^b-i3) = np + np - - V??- (A24) 

Combining (A22) and ( |A24| |, noting that the remainder terms ±(nB + np) cancel one-another, we obtain 

(fiB + nF){N + fiB — np) — + 29^9 — 2fi^ri , (A25) 


so we can identify; 


where 


=2Q{N+1-Y), 


Q = hp + fiB, 

' ^ ^ — 29^9 + 27 )^ 7 ] 

Y = hF - UB + 1 + ^- —■ 

Q 


(A26) 

(A27) 

(A28) 


5. Sum rule for spin and super-Casimir 

The second Casimir invariant of fho SU{2\1) group can be writterP^ 

X^=TT[XmX], (A29) 

where X = Xab is the three-dimensional matrix formed out of the Hubbard operators and m = diag(l, 1,-1). If we 
expand this result we obtain 


= XalsX^a — [^aO.^Oa] — (-^OO)^, (A30) 

with an implied summation over the repeated indices [3 = ±. Substituting for the Hubbard operators using ( p!^ 
we obtain 


x' = + + a!_ - X^o + {A+_, A_+} - [A+o, Ao+] - [A_o, Ao_] 

= [{tip - 77b)/2]'^ + [(A -Tip- 77b)/ 2f - flp + '*/’} “ ^]/2 - [??,f?V2 

= A^/4 - (tip + nB){N - Tip + 77 b)/ 2 + {7p', 7p} - [9\ 9]/2 + [nf ,7f\/2. 


(A31) 
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Using the Hubbard algebra of the SU{2\1) generators we have 

= 2V’V + bP, V’1 = + [^-+, ^+-] = 2V’V + ^— - ^++, 

[0\ 0\I2 = 0^0 - {0, 0t}/2 = 0^0 - {Xo+, X+o} = 0^0 - (Xoo + X++), 

[r]\ 7?]/2 = r}<r] - {?7, r]^} = rj^r] - {X_o, Xq-} = rj^ri - (X__ + Xqo), 

Adding up these expressions, we find that 


{ip^, ?/)} — [6*'^, 0]/2 + [t]^, ri\/2 = 2'tp^xp — 0^0 + rj^r]. 


Inserting this into (A31) we obtain 

= N"^/A. — {np + nB){N — np + ub )/2 + 2xp'''xp — 0^0 + 


Finally, comparing with {A25), we obtain 


= N^A-S^/A. 


The corresponding sum rule 


N^/4 = x^ + S^/A. 


(A32) 

(A33) 

(A34) 

(A35) 

(A36) 


expresses the fact that sum of the pair/charge fluctuations and the spin fluctuations is a constant. 


Appendix B: Computation of the Fermionic part of 
the free energy 

The fermionic part of the solution reduces to two de¬ 
coupled impurity problems. Here we explicit the calcu¬ 
lation for a single impurity. The partition function for a 
single impurity can be written as: 

Zf = J VfiFe~^^, Viip =V[c, f.v.Xp], 

already transforming from imaginary time to Matsubara 
frequencies, the action can be written as: 


is the f-fermion free energy, where 

is the bare conduction electron propagator. We evaluate 
the sum over k in Sj as an integral over energy with a 
constant density of states. Analytically continuing the 
Matsubara frequencies to the real axis oo ^i5)\ 

k 

= kl^y p(e)de ^ ^ =F iTT5{w - e)^ , 

= -iTe{D - \uj\)sgn{uj), (B4) 


Sp= ^ ^ + <^k)Ck<T 

n L kcr 

<j 

+ 'Z + h.C. j 


(Bl) 


We start by integrating out the conduction electrons, 
taking into account their effect in the self energy of the 
fermions that compose the spin. The effective fermion 
propagator can be written as: 

G/ = [(G°/)-'-%]-\ (B2) 

where — Xf is the bare f-fermion propaga¬ 

tor, and 


where T = Trpol'^P, D the bandwidth, po the constant 
density of states, &{x) the Heaviside step function. Here 
uj indicates the imaginary part of the frequency. 

Now the fermionic part of the free energy reads: 

Sf — ^ ^ /or(^^n)( T ^F T f 

na 

AluP 

- PNXfQf, (B5) 

Jk 

where T^ = TQ{D— \ujn\)sgn{ujri)- Now we can integrate 
out the f-fermions and write an effective action at the 
saddle point values of v and Xf (to be determined by 
extremization of the free energy, see main text): 

S^ff ^ ^ log[-iw„ + Xf + iTn] 

na 

- pNXpqp. (B6) 
Jk 


k 


(B3) 


The sum over Matsubara frequencies can be performed 
as an integral in the complex plane weighted by the Fermi 
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distribution function f{z) = — 1) Note that the 

integral involves a branch-cut: 

^ log[-iuin + Xf + iTO{D - \uin\)sgn{Un)] 

na 

^ ^ Jc + Af + ire{D - \z\)sgn{z}]f{z) 


2711 


+ 


f dzf{z) log[-z + Af - ir] 
J-D 

[ dzf{z) log[- 2 : + Af + ir] 
Jd 


(B7) 


which simplifies to: 

rD 


— f dzf{z)Im[\og[—z + XF + iT]]‘ (B8) 

^ J-D 

In the zero temperature limit the Fermi function sets 
the upper limit of the integral to zero. Evaluating the 
integral we find the free energy: 


Ff 1 
^ = -Im 
N TT 


{Xf + 'i^F) In 


Ai? “h iF 

De 


- - XfQf, 

Jk 


that can be rewritten as: 


Ff 1 
— = -Im 
N IT 




once we define 

= Ai? + iF, 

and the Kondo temperature 

Tk = . 


(B9) 

(BIO) 

(Bll) 

(B12) 


Appendix C: Computation of the bosonic part of the 
free energy 


) . (C3) 

Integrating out the bosons and taking the saddle point 
value of Ab and g, which will be determined by the ex- 
tremization of the free energy with respect to these pa¬ 
rameters, we can write: 


oEff 

Zb = e B ^ 


(C4) 


where 


= F.^og[Det[LBiiiyn)]]+^N^-j^ 


-2/3NXBiqB + l/2), (C5) 

n,x=± 

—2/3NXB{qB + 1/2), 


Jh 


where 


E% = ±JXl-\g\F 


(C6) 


The sum over Matsubara frequencies can be written in 
terms of an integral over the imaginary plane weighted 
by the bosonic distribution function n{z) = — 1)“^: 


n,x=± 


dz 


^ \og[E%-mn] =j^^,\og[E%-z\n{z). 

(C7) 

In the zero temperature limit: 

Y^og[Det[LB{iVn)]] ^ NYi-E%M-E%), 

(C8) 


X = ± 


where &{x) is the Heaviside step function, so that the 
bosonic part of the free energy reads: 

^ = \/>^FF+Z-2XB{qB + l/2). (C9) 


From the main text we have that the bosonic part of 
the partition function is: 

Zb = J 'DiiBe ~^^, V^b = 9, Xb] , 

already transforming from imaginary time to Matsubara 
frequencies, the action can be written as: 

Sb = Y^B<Ti'^^n)LB{il'n)'^Ba{il'n) 

no- 

+/37vkl! _ 2/37 VAb(</b + 1/2), (Cl) 

Jh 


where 




Appendix D: Computation of the fluctuations of the 
fermionic hybridization 

In this appendix we define and compute Xcb{^^n)‘ 
From the main text we have: 

Xcbij^r) — ~n ^ ^ (Bl) 


kk'm 


where, 


GbiiVn) = {iVn - is) ^ (D2) 

is the bosonic propagator, with ^b = \/Ag — |(/p, and 

G\z\z' {iuJn) ~ G^(^iuJn)dk.k.' 

+ \v\‘^G^{iuJn)Gf{iuJn)Gl,{iuJn), (D3) 
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with the propagators defined in Appendix [B| 

Evaluating the sum over momenta: 

{iLdn) = 

kk' k 

+ 5 ( 04 ) 

where 

^ Gk(^^n) = -i7rposgn{uJn), (D5) 

k 

as computed in the evaluation of the fermionic part of 
the free energy. In the infinite bandwidth limit the sum 
over k can be written as: 


y^Gkk^ {iuJn) = -i7Tposgn{uJri) 

kk' 

_ ttPoE _ 

iujn — Xf iTsgn{uJn) 


(D6) 


where po is a constant density of states and E = 7rpo\v\‘^ 
as before. 


Back to the computation of Xcb- 

Xcb{i^r) = Xlbii^^r) + Xlbii^^r), 

where 

1 1 ^ mposgn{wm) 

2fi J 

_ iTTpo 
2Fi 

rD 


+ 


Z -\-iuJr - 
[27ri/(AB - iu)r)i-l) 
(- 1 ) 


[ dzf{z) . 

J^D Z + lWr- is 

/■■”*/(,)_ (+1) 

JD 


Z + lLdr - 


(D7) 


(D8) 


where z = Im{z) and /{Xb —ioJr) = —ti{Xb)- In the zero 
temperature limit f{z) 0(—z) and n{XB > 0) ^ 0, so: 

The second part of Xcbi'^^r)- 


Xcbii^r) = ^^-— 


KpoT 


p ^ ibJm - \f+ iTsgn{u)m) + iojr - is 


(DIO) 


can be computed in analogous fashion; 


Xlbii^r) = 


por 


2i iujj. — ^B T Xf T 
Poh _1_ 

2i iujr ~ ^B Xf — iT 


Log 

Log 


-Xf ~ iT 


— Xf — — D 

—Xf + 

—A + ir — D 


- Log 
- Log 


-^B + 


-^B +iuJr - D 

—^B + 

—^B + 'i'^r ~L>J 


(Dll) 

Continuing to real frequencies ujr ^ UJ — iS, writing 1/Jr = —poLog\{X + iT)/D\, in the infinite bandwidth limit: 


Xcb{^ -iS) - ^ = ^poLog 
Jk 

PqL 


Ai? + ir 


£,b - ^ 

1 


iFpQQ{uj — A) + 


Poh 


1 


2i uj — i5 — ^B Xf — iL 


^Log 


2i uj — iS — ^B + Ai? + iT 
A^ — ir 


^Log 


Ai? + ir 

^B - bJ -^iS 


^ 


(D12) 


In the transition line, where ^b = ^ V^% ~ 9^ = Xb = Xf = X^ we have the simplified form: 


Xcb{oj -i5) - ^ = -LpoLog 
Jk 


A Air 


A — UJ i5 

por 1 


2i cj2+r2 


— iFp{)Q{uj — X) 

A + ir 


{uj — iV)Log 


X — UJ i5 


— (cj + iV)Log 


X-iV 
X — UJ i5 


rewriting, 


XcbUj -iS) - — = potoRe 
Jk 


Log{X + iV) 


UJ T iT 


Po^ 

r2+a;2 


log(A - UJ + i6), 


(D13) 

(D14) 
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which is the form of Xcb discussed in the main text. 
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